本发明涉及一种基于对由碰撞产生的声信号的分析来确定物体表面上的碰撞位置的方法。该方法还包括信号处理步骤,在所述信号处理步骤中,对声信号进行加权,以具体由于在物体边界处的反射而引起的假性贡献。
1.一种用于确定物体(1)表面(9)上碰撞位置(11)的方法,基于对由碰撞产生的声信号的分析来确定物体表面上碰撞位置,其特征在于:还包括信号处理步骤,在所述信号处理步骤中,对声信号进行加权,从而考虑具体由于物体(1)边界处的反射而引起的假性贡献,其中信号处理步骤包括识别声信号的第一波前的步骤以及应用阻尼窗的步骤,其中应用阻尼窗的步骤包括使用声信号与以下指数阻尼窗相乘:ntent="drawing" img-format="tif" inline="no" orientation="portrait" wi="700"/>其中t_i与第一波前的到达时间相对应,Δt是时间延迟,以及T是指数时间常数。
2.根据权利要求1所述的方法,其中在时域执行所述加权。
3.根据权利要求1所述的方法,其中识别第一波前的步骤包括识别最大信号幅度。
4.根据权利要求3所述的方法,其中识别步骤还包括:将第一波前的到达确定为当声信号超过基于最大信号幅度的阈值时的时间。
5.根据权利要求1至4中任一项所述的方法,其中由至少两个变换装置(5a-5d)来感测声信号,针对一个变换装置,来识别第一波前使得针对所有的变换装置应用阻尼窗的相同时间原点。
6.根据权利要求5所述的方法,其中只针对一个变换装置来识别第一波前,使得针对所有的变换装置应用阻尼窗的相同时间原点。
7.根据权利要求1至4中任一项所述的方法,其中由至少两个变换装置(5a-5d)来感测声信号,使用由所述至少两个变换装置感测到的信号的互相关结果来执行声信号的加权。
8.根据权利要求7所述的方法,其中识别第一波前贡献的步骤包括确定所述互相关结果的最大值。
9.根据权利要求1所述的方法,其中应用阻尼窗的步骤包括应用对称窗函数。
10.根据权利要求9所述的方法,其中应用阻尼窗的步骤包括应用高斯窗、汉明窗、汉宁窗或布莱克曼窗。
11.根据权利要求1至4和9-10中任一项所述的方法,还包括:对声信号应用低通滤波和/或高通滤波。
12.一种用于基于确定物体表面(9)上碰撞位置(11)的装置,基于对由碰撞产生的声信号的分析来确定物体表面上碰撞位置,所述装置包括:信号处理装置,对声信号进行加权,从而考虑具体由于物体(1)边界处的反射而引起的假性贡献,其中信号处理装置包括识别声信号的第一波前的装置以及应用阻尼窗的装置,其中应用阻尼窗的装置包括使用声信号与以下指数阻尼窗相乘的装置:ntent="drawing" img-format="tif" inline="no" orientation="portrait" wi="700"/>其中t_i与第一波前的到达时间相对应,Δt是时间延迟,以及T是指数时间常数。
13.根据权利要求12所述的装置,其中在时域执行所述加权。
14.根据权利要求12所述的装置,其中识别第一波前的装置包括识别最大信号幅度的装置。
15.根据权利要求14所述的装置,其中识别装置还包括:将第一波前的到达确定为当声信号超过基于最大信号幅度的阈值时的时间的装置。
16.根据权利要求12至15中任一项所述的装置,其中由至少两个变换装置(5a-5d)来感测声信号,针对一个变换装置,来识别第一波前,使得针对所有的变换装置应用阻尼窗的相同时间原点。
17.根据权利要求16所述的装置,其中只针对一个变换装置来识别第一波前,使得针对所有的变换装置应用阻尼窗的相同时间原点。
18.根据权利要求12至15中任一项所述的装置,其中由至少两个变换装置(5a-5d)来感测声信号,使用由所述至少两个变换装置感测到的信号的互相关结果来执行声信号的加权。
19.根据权利要求18所述的装置,其中识别第一波前贡献的装置包括确定所述互相关结果的最大值的装置。
20.根据权利要求12所述的装置,其中应用阻尼窗的装置包括应用对称窗函数的装置。
21.根据权利要求20所述的装置,其中应用阻尼窗的装置包括应用高斯窗、汉明窗、汉宁窗或布莱克曼窗的装置。
22.根据权利要求12所述的装置,还包括:对声信号应用低通滤波和/或高通滤波的装置。技术领域
本发明涉及一种基于碰撞产生的声信号的分析来确定物体表面上碰撞位置的方法。
背景技术
根据WO2006/108443已知这种方法,WO2006/108443公开了一种确定包括N个换能器的表面上的碰撞位置的方法,将感测信号S_i(t)传输给处理单元,其中所述方法包括以下步骤:a)计算所述N个不同换能器的感测信号的傅里叶变换的P个互相关结果(inter-correlationproducts);b)计算P个互相关结果的傅里叶逆变换p’_ij(u);c)针对所述物体表面的每一个区域k计算P_k(u)=∑p′_ij(u-τ_ijk);以及d)找到区域k₀,其中,在所有区域k的特征值P_k(u)中,区域k₀的特征值P_k0(u)最大。这里,τ_ijk对应于由每一个区域k的两个不同换能器感测的信号的到达时间差。因此,定位方法基于在发生碰撞位置处交叉相关函数(cross-correlationxa0function)和的显著峰的存在。然而在特定情况下,由于物体的较宽区域中的声信号反射或者物体中的其他不连续性(例如螺旋、裂缝或孔洞)而引起的假性贡献使得碰撞位置的确定出错。已经提出了通过在预定时间间隔之后切断信号或者通过使用低通滤波器来解决这一问题。然而,所提出的方案获得的结果并不令人满意。在接收到的没有经过反射的感测信号与反射后的感测信号之间不可能清楚分离的情况下(例如这可能发生在由于材料性质导致信号分散的情况下),所述切断导致修正的信号,由此也限制了碰撞的定位精确度。
发明内容
因此,本发明的目的是提出了一种考虑到假性贡献(例如物体边界处反射的贡献)来更好地处理感测信号的方法。利用根据权利要求1所述的方法实现了该目的。通过对声信号进行加权,可以保持所述信号的完整信息,同时,相对于由物体边界处的反射引起的贡献(典型地随后到达),放大从碰撞位置直接接收的感测信号的贡献。因此,碰撞位置的确定变得更加可靠和精确。在上下文中,应该理解,所处理的声信号不必须是由换能器立即感测的信号,而可以是已经过例如模拟数字变换处理的信号。优选地,在时域执行加权,由此在识别碰撞位置方面实现最佳的结果。根据优选实施例,信号处理步骤可以包括:识别声信号的第一波前的步骤;以及应用阻尼窗的步骤。碰撞之后声信号在物体上展开,并且由特定数量的变换装置(例如换能器)来感测所述声信号。因为由于几何原因,只有在声信号的第一波前已经到达正讨论的换能器之后,才能感测信号的第一反射贡献,所以第一波前的识别是鲁棒参数,以识别信号的一部分,该部分应该具有比感测信号的其余部分的加权因子高的加权因子,并且应该被应用阻尼窗的步骤所考虑。在上下文中,术语“第一波前”意味着在没有因反射而造成假性贡献的情况下到达变换装置的信号部分。优选地,识别第一波前的步骤可以包括识别最大信号幅度。所述最大信号幅度是易于检测的参数,因此不要求很大的计算能力,同时是用于找出声信号开始的鲁棒的起始参数。感测声信号的备选参数可以是信号的功率、最大峰到峰幅度或者信号的能量,而不是最大幅度。在识别信号幅度之前,可以执行信号整流和平滑的另外步骤以促进数据分析。根据有利的实施例,识别第一波前的步骤还包括:将第一波前的到达确定为当声信号超过基于最大信号幅度的阈值时的时间。例如,可以将声信号的第一波前的到达归结于达到最大幅度的特定百分比的时间。优选地,应用阻尼窗的步骤可以包括使用声信号与以下指数阻尼窗的相乘:ntent="drawing" img-format="tif" inline="no" orientation="portrait" wi="700"/>其中t_i与第一波前的到达时间相对应,Δt是时间延迟,以及T是指数时间常数。这看来针对这种指数阻尼窗在数据分析中达到最佳结果,其中项Δt(Δt>0)是将阻尼窗的起始优选地移动至较早时间以优化结果的时间延迟。根据应用,Δt可以是固定值或自适应的,例如根据装置的初值而自适应。有利地,可以由至少两个变换装置感测所述声信号,并且其中针对一个变换装置,具体地只针对一个变换装置,来识别第一波前,使得对所有的变换装置应用阻尼窗的相同时间原点。由此,可以保持计算量较低。使用第一感测信号的第一波前获得最佳结果。在上下文中,通过对N个信道的信号幅度使用固定阈值来识别第一感测信号。根据优选实施例,可以由至少两个变换装置来感测声信号,并且可以使用由所述至少两个变换装置感测的信号的互相关结果来执行声信号的加权。因此与上述不同,没有对直接感测的信号进行加权,而是对互相关进行加权,优选地将所述互相关用于确定碰撞的位置。再次地,如上述方法一样,可以在不会丢失来自直接信号的信息的的情况下抑制反射的不利影响。为了更进一步地减小假性贡献的影响,本发明的变体包括:如上述公开的在第一步骤中对感测信号进行加权,以及在第二步骤中对交叉相关(cross-correlation)结果进行加权。有利地,识别第一波前的步骤可以包括确定互相关结果的最大值。作为备选,也可以使分析基于对于最大平方幅度、最大峰到峰幅度、均方根等等。这些是用于识别两个换能器之间的信号到达时间差作为第一级近似的鲁棒参数,这两个换能器的互相关已确定。优选地,应用阻尼窗的步骤可以包括应用对称窗函数,具体地应用高斯窗、汉明(Hamming)窗、汉宁(hanning)窗或布莱克曼窗(Blackman)窗,以提高直接信号的贡献,以及在不会完全丢失内容的情况下使反射的贡献衰减。根据有利的变体,所述方法还可以包括向声信号应用低通滤波和/或高通滤波,这导致可以从中得出碰撞位置的数据的进一步优化。所述方法还涉及一种包括一个或多个计算机可读介质的计算机程序产品,所述计算机可读介质具有计算机可执行指令,所述计算机可执行指令用于执行任一前述权利要求中所述方法的步骤。还利用根据权利要求13的装置来实现本发明的目的,利用该装置可以达到与所述方法相同的有利效果。
附图说明
以下将结合附图描述本发明的有利实施例,其中:图1是使用根据本发明的确定碰撞位置的方法的触敏装置的示意性三维视图;图2示出了用于确定碰撞位置的交叉相关结果;图3示出了当k与k0不同时特征值P_k(u)较低的事实,其中k₀是碰撞开始的区域;图4示出了变换装置感测到的声信号的反射的效果;图5示出了反射对于确定碰撞位置的精确度的不利影响;图6示出了感测信号,所述感测信号已经被处理以识别声信号的第一波前;图7示出了用于对感测声信号进行加权的阻尼窗;图8示出了用于对交叉相关结果进行加权的阻尼窗;以及图9示出了加权结果。
具体实施方式
图1是使用根据本发明的确定碰撞位置的方法的触敏装置1的示意性三维视图。触敏装置1包括:交互装置3、变换装置5a-5d(变换装置5d在图1中不可见)和与所述变换装置5a-5d链接的信号处理装置7。交互装置3是透明的、半透明的或者不透明的面板,所述面板是平坦的或者是弯曲的,并且可以由诸如皮革、乳胶、硅、塑料、玻璃、金属、木头、石膏板或复合材料之类的任意合适材料制造。它也可以是刚性的或者是柔软的,只要声信号可以在其中或其上传播。交互装置3提供用于允许在用户和装置1之间进行交互的表面。在该实施例中,用户典型地在与变换装置相反的表面侧9上提供碰撞,而变换装置5a-5d也可以被布置在发生碰撞的那一侧上。在图1中,示出了具有四个变换装置5a-5d的触敏装置1。然而,变换装置5a-5d的数目并没有限制,只要至少两个变换装置作为所述触摸控制系统1的一部分,就可以实现本发明的有利效果。变换装置5a-5d包括换能器,如,压电换能器形式的应变仪,所述压电换能器例如是陶瓷膜或PVDF膜。然而,变换装置5a-5d不局限于应变仪的使用,因为通常被动的或被供电的任意类型的压力或力传感器都是合适的,只要它可以将碰撞产生的声信号转变成电信号。变换装置5a-5d将感测信号发送给信号处理装置7,所述信号处理装置被配置为分析感测信号,使得可以确定碰撞位置11。另外,处理装置7可以被配置为向附加的装置输出碰撞位置11,或者触发与物体上发生碰撞的位置相链接的动作。信号处理装置7包括信号调节单元,用于例如通过对来自变换装置5a-5d的感测信号进行滤波、放大和/或转换为数字信号来适配该感测信号。该信号然后可以由碰撞定位单元进一步处理,所述碰撞定位单元处理由交互板3上的碰撞引起的数字信号,以确定交互板3上的该碰撞的位置(例如x、y坐标)。现在将描述确定碰撞位置的一种方法:首先,在碰撞之后,四个声学换能器5a-5d提供感测信号s_i(t),i从1到4,然后信号处理装置7(信号调节单元)将感测信号转变为离散的感测信号s_i(n)。之后,执行傅里叶变换以获得针对每一个形成装置5a-5d的S_i(w)。信号调节单元最终也考虑分散效应。然后,对于所有可能的组合执行交叉相关计算:P_ij(w)=S_i(w)S_j(w)*。使用傅里叶逆变换将所述交叉相关重新变换到时域。在下一个步骤中,对各交叉相关结果求和p_k(u)=∑_ijp_ij(u-τ_ijk),其中τ_ijk对应于针对交互板3表面上给定位置k的到达时间差。到达时间差值τ_ijk与由位置k处的碰撞产生的声信号相对于装置5a-5d中变换装置i和j的位置的传播时间差相对应。所述值是理论上确定的,或者处于装置1的起始点或校准点。然后为了识别碰撞的位置11,查找p_k(u)的能量最大值。代替能量,也可以使用交叉相关和的其他参数,例如最大值、功率或最大平方幅度、最大峰到峰幅度、均方根等等。可以在上述现有技术文献中找到有关该方法的另外细节。图2示出了所述方法的理论结果。在左侧,针对给定个数的换能器示出了交叉相关结果,在中间,通过针对与碰撞位置11相对应的位置k₀的到达时间差值τ_ijk对所述交叉相关进行校正,右侧示出了针对位置k₀的具有预期显著峰值的交叉相关结果之和。图3现在示出了当k不与k₀相对应时的特征值p_k(u)。在这种情况下,所述特征值p_k(u)远小于p_k0(u)。然而实际上,对于感测信号的假性贡献的存在(例如在交互片3的边界处反射的存在)使情况恶化。图4的上部分示出了在自由场传播的理论情况下的感测信号的一个示例。在没有发射的情况下,信号的持续时间在时间上有限。现在在由于非理想边界条件导致存在反射的情况下,信号在时间上展开(图4的下部分)。如图4的下半部分所示,不再可能区分在边界上没有反射的情况下直接到达变换装置的第一波前和由反射引起的第一贡献。事实上,由于交互板上的碰撞,声信号能够已经具有一定的持续时间,使得直接感测的信号的末端可以干扰反射信号的起始。另外,分散效应进一步模糊了所述信号。在包括假性贡献的感测信号上执行上述确定碰撞位置的方法,获得了如图5所示的结果。三维曲线示出了p_k(u)作为交互板3上的位置k的坐标x、y的函数。由于虚假信号,与碰撞位置11相对应的交叉相关结果之和51的最大值不如图2所示的显著。此外,假性贡献错误地加强了不与实际碰撞位置k₀相对应的一些其他位置53、55、57、59的交叉相关结果。因此,更加难以正确地识别碰撞的位置11。在最差情况下,信号处理装置7甚至可能识别错误的碰撞位置。为了克服这一问题,根据本发明的方法提出执行感测声信号的加权,使得在不完全去除信号尾部的情况下,与末端信号相比,起始点处的贡献具有更大影响。在模拟域(因此,在执行定位例程之前,在信号调节单元内)中或者在数字域中,在信号处理装置7中执行数据处理。本发明的信号处理步骤包括:第一步骤,识别声信号的第一波前;以及第二步骤,应用阻尼窗。第一波前的识别基于在每一个变换装置5a-5d处的信号到达时间的估计。该实施例中到达时间估计是通过以下步骤来实现的:首先,执行全波整流;然后例如通过对完全整流的信号执行一阶低通滤波来执行平滑例程。一旦已经估计了信号包络,针对每一个变换装置5a-5d确定了最大幅度a_i^max。所述最大值不必须与第一波前的到达时间相对应,已经发现所述最大值与可以从中估计到达时间的鲁棒参数相对应。为此目的,确定信号何时第一次达到该最大级别的给定百分比。然后认为该时间点表示声信号到达变换装置的到达时间。因此,所述到达时间与信号达到α×a_i^max的时间n_j(按照离散信号形式(n))相对应,其中0<α<1。当α在0.15至0.20的范围内(优选地是0.15)时,实现了最佳结果。图6示出了所述情况,示出了完全整流且平滑后的信号,其中也指示了最大幅度A_i^max和与α×A_i^max相对应的时间ni处的到达时间。应用阻尼窗的第二步骤包括将感测信号与合适的阻尼窗相乘,优选地与如图7虚线所示的指数窗相乘。所述指数窗可以以与到达时间相对应的时间索引n_i开始,然而当阻尼窗比到达时间更早开始时实现了最佳结果。因此在指数窗的情况下,可以将阻尼窗定义为:ntent="drawing" img-format="tif" inline="no" orientation="portrait" wi="700"/>这里Δt与正值相对应,T是指数时间常数。对于离散时间索引,等式变成:w(n)=exp((n_i-Δn-n)/N)。其中离散指数时间常数N=T/Ts,Ts是采样时期,T是指数时间常数。从图7可以看出,通过应用阻尼窗,不是仅在特定持续时间之后切断信号,而是在整个信号持续时间上保留信息。根据该实施例的变体,可以将到达时间的确定只局限于变换装置之一,具体地局限于首先观察到信号的变换装置。由此,可以限制计算量,并且阻尼窗可以对于所有感测信号在相同的时间t_i或者t_i-Δt开始。根据另外的变体,加权公式也可以是:ntent="drawing" img-format="tif" inline="no" orientation="portrait" wi="700"/>在这种情况下,没有取消感测信号的开始。指数阻尼窗只是一个发明性示例,也可以使用其他类型的阻尼窗。实际上,可以使用任意类型的阻尼窗,阻尼窗可以是非对称的,例如线性的,或者可以是对称的,例如高斯、布莱克曼、汉宁或汉明。在对称窗的情况下,可以以t_i-Δt附近为中心。这可以防止无意地消除感测信号上的第一波前。根据本发明的第二实施例,在变换装置5a-5d感测到的信号的互相关结果上执行声信号加权方法。在本发明方法的第二实施例中,因此可以识别互相关结果的最大值,阻尼窗用于提高由直接信号引起的交叉相关的贡献,以及用于使反射产生的信号的贡献衰减,优选地所述阻尼窗是对称窗,例如高斯、汉明、汉宁或布莱克曼窗。图8示出了该第二实施例的概念。示出了由两个感测信号(平纹(plain)线)计算出的交叉相关结果。虚线的阻尼窗以交叉相关结果的最大值为中心。另外,交叉相关的最大值提供了两个不同变换装置的信号飞行时间差的良好一阶近似。这里高斯窗如下定义: w ( n ) = e - 1 2 ( n - ( N - 1 ) / 2 σ ( N - 1 ) / 2 ) 2 , σ ≤ 0.5 ]]>其中,在44kHz的采样频率下,σ=1e-2并且N=2048。图9说明了用于识别碰撞位置的交叉相关结果的加权的效果。可以对本发明的第一和第二实施例进行组合以进一步抑制假性贡献。还可以通过应用高通滤波器和/或低通滤波器去除例如由电源等引起的对信号的假性贡献,来进一步改进第一和第二实施例。利用根据本发明的第一和第二实施例,可以较大程度地抑制例如由交互装置边界处的反射引起的假性贡献的不利影响,而不会丢失确定碰撞位置所必须的信息。